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ABSTRACT 


Recent  conflicts  have  highlighted  the  difficulties  in  using  aircraft  to  supply  troops 
in  modem  warfare.  Lift  fan  technology  is  seen  as  one  way  in  which  to  improve  future 
supply  aircraft  to  meet  the  needs  of  the  military.  A  previous  study  designed  a  heavy  lift 
aircraft  with  lift  fan  engines  that  used  future  engine  technology.  This  present  study 
modified  the  design  by  using  current  engine  technology  for  the  lift  fan  engines.  The 
modification  is  important  because  a  design  that  uses  current  technology  is  more  likely  to 
be  brought  into  service  in  the  near  future. 

This  thesis  documents  the  process  required  to  use  current  technology  in  a  lift  fan 
engine  for  a  heavy  lift  aircraft.  The  process  uses  current  software  and  focuses  on  the 
design  of  the  following  components:  the  powerplant,  the  transmission  shafts  and  the  lift 
fan.  The  result  is  a  propulsion  system  which  allows  a  185,000  lb  aircraft  to  takeoff 
vertically  as  well  as  cmise  at  speeds  greater  than  Mach=0.6. 


V 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


VI 


TABLE  OF  CONTENTS 


1.  INTRODUCTION . 1 

A,  BACKGROUND . 1 

B,  COMPONENTS . 1 

C,  BASELINE  AIRCRAFT . 1 

D,  PERFORMANCE  ESTIMATES . 2 

1,  Background . 2 

2,  Results . 2 

E,  ENGINE  PERFORMANCE  CALCULATIONS . 3 

F,  INSTALLED  ENGINE  PERFORMANCE . 4 

IL  POWERPLANT . 5 

A.  DESIGN  CONSIDERATIONS . 5 

B.  BASELINE  ENGINE . 6 

C.  CHANGES  TO  THE  BASELINE  ENGINE . 6 

D.  DERIVATIVE  ENGINE . 7 

1,  Referencing  Engine  Performance  to  New  Design  Point . 7 

2,  Turbine  and  Nozzle  Changes . 8 

3,  Engine  Changes  with  Changes  in  LP  Mechanical  Efficiency . 9 

4,  Selection  of  the  Derivative  Engine  Design  Point . 12 

E.  ENGINE  PERFORMANCE  IN  CRUISE . 13 

1,  Design  Considerations . 13 

2,  Cruise  Performance . 14 

F.  BLEED  AIR . 16 

III.  TRANSMISSION  SHAFTS . 19 

A.  DESIGN  CONSIDERATIONS . 19 

B.  FAN  AND  SHAFT  ROTATIONAL  SPEED . 19 

C.  TRANSMISSION  SHAFT  THICKNESS . 20 

1.  Governing  Equations . 20 

2.  Shaft  Design . 21 

C.  TRANSMISSION  SHAFT  LENGTH . 22 

1.  Critical  Speed . 22 

2.  Governing  Equations . 22 

3.  Actual  Design . 23 

IV.  LIFT  FAN . 25 

A.  DESIGN  CONSIDERATIONS . 25 

B.  AEROTHERMODYNAMIC  ANALYSIS . 25 

C.  LIFT  FAN  PRELIMINARY  DESIGN . 28 

1.  Design  Tool . 28 

2.  Determination  of  Initial  Lift  Fan  Paramaters . 29 

a.  Fan  Blade  Tip  Speed  (Ut) . 29 

b.  Inlet  Relative  Flow  Angle  (fii) . 29 

vii 


3.  Design  of  Lift  Fan . 30 

a.  Design  Iteration . 30 

b.  Lift  Fan  Gearing . 30 

c.  Fan  Designs. . 31 

V.  FINAL  AIRCRAFT  DESIGN . 35 

A.  ADDITION  OF  THRUST . 35 

1.  Requirement . 35 

2,  Solution . 35 

B.  DESIGN  OF  TRANSMISSION  SHAFT . 36 

C.  FINAL  AIRCRAFT  DESIGN . 37 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS . 39 

A,  CONCLUSIONS . 39 

B,  RECOMMENDATIONS . 39 

LIST  OF  REFERENCES . 41 

APPENDIX  A . 43 

APPENDIX  B . 47 

COMPONENT  MAPS  FOR  0  FEET  AND  M=0.0  (CRUISE) . 47 

COMPONENT  MAPS  FOR  5000  FEET  AND  M=0.55  (CRUISE) . 48 

COMPONENT  MAPS  FOR  20000  FEET  AND  M=0.6  (CRUISE) . 49 

APPENDIX  C . 51 

INITIAL  DISTRIBUTION  LIST . 59 


viii 


LIST  OF  FIGURES 


Figure  1 .  Two  Spool  Mixed  Flow  Turbofan  Engine  Configuration . 4 

Figure  2.  Variation  of  Engine  Thrust  with  A8 . 11 

Eigure  3.  Variation  of  Total  LP  SHP  with  A8 . 11 

Eigure  4.  GASTURB  Eimiters  Seleetion  Sereen . 13 

Eigure  5.  GASTURB  Modifiers  Selection  Screen . 14 

E igure  6 .  Engine  Ean  Map . 15 

Eigure  7.  Engine  EP  Turbine  Map . 16 

Eigure  8.  Plot  of  Shaft  Eength  Versus  Natural  Erequency . 24 

Eigure  9.  Mass  Elow  vs  Thrust . 27 

Eigure  10.  Mass  Elow  vs  Eift  Ean  Diameter . 28 

Eigure  11.  Basic  Blade  Theory . 29 

Eigure  12.  Eift  Ean  Engine  Installation  Dimensions . 36 

Eigure  13.  Top  View  of  Heavy  Eift  Aircraft . 37 

Eigure  14.  Perspective  View  of  Heavy  Eift  Aircraft . 38 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


X 


LIST  OF  TABLES 


Table  1.  Aircraft  Characteristics . 2 

Table  2.  Aircraft  Thrust  and  Fuel  Requirements  (Aircraft  Weight:  180000  lbs) . 3 

Table  3.  Estimated  Baseline  Engine  Parameters . 6 

Table  4.  Estimated  Eift  Engine  Parameters . 8 

Table  5 .  Variation  of  Engine  Parameters  with  EP  Mechanical  Efficiency . 10 

Table  6.  Estimated  Eift  Engine  Parameters . 12 

Table  7.  Cruise  Performance  for  One  Engine . 15 

Table  8.  Estimated  Derivative  Engine  Parameters . 17 

Table  9.  Properties  of  Titanium  Alloy  (Ti-6A1-4V) . 21 

Table  10.  Values  Used  for  Calculation  of  Maximum  Shear  Stress . 21 

Table  1 1 .  Values  for  an  for  Clamped-Clamped  Beam . 23 

Table  12.  initial  Eift  Fan  Parameters . 27 

Table  13.  Eift  Ean  Design  Constraints . 30 

Table  14.  Diameter  vs  Rotational  Speed  (Ean  Blade  Tip  Speed=1632ft/s) . 3 1 

Table  15.  Potential  Eift  Ean  Designs . 32 

Table  16.  Detailed  Ean  Design  Parameters . 33 

Table  17.  Summary  of  Thrust  Calculations . 35 

Table  18.  P-35C  Engine  Performance . 35 

Table  19.  Summary  of  Thrust  Calculations . 36 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


ACKNOWLEDGMENTS 


I  would  like  to  thank  the  following  individuals  for  their  guidanee  with  this  thesis: 


Dr.  Raymond  Shreeve,  Professor,  Naval  Postgraduate  School 
Dr.  E.  R.  Wood,  Professor,  Naval  Postgraduate  School 

Dr.  Paul  Bevilaqua,  Manager  of  the  Advanced  Developmental  Programs  Division, 
Lockheed  Martin,  Palmdale,  CA 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


XIV 


I.  INTRODUCTION 


A,  BACKGROUND 

Recent  conflicts  have  highlighted  the  difficulties  of  supplying  troops  in  modern 
warfare.  A  problem  with  staging  supplies  in  foreign  countries  has  led  the  military  to 
examine  supplying  troops  directly  from  ships.  A  recent  study  at  the  Naval  Postgraduate 
school  outlined  the  requirements  that  a  future  cargo  aircraft  might  face  in  this  scenario: 
large  enough  to  carry  37,500  pounds  of  cargo  including  vehicles;  be  able  to  cruise  at 
relatively  high  speeds  on  missions  up  to  600  nautical  miles;  and  be  flexible  enough  to 
land  and  takeoff  vertically  (Erhardt,  W.,  2002).  One  proposal  for  this  mission  involved  an 
aircraft  that  used  lift  fan  engines  (Wood,  E.  R.,  2003).  The  present  project  attempted  to 
improve  on  the  aircraft  outlined  in  the  proposal  by  incorporating  lift  fan  engines  that 
require  only  modest  changes  to  the  engine  designs.  The  approach  taken  is  outlined  in  the 
following  paragraphs. 


B,  COMPONENTS 

The  present  project  addressed  the  following  components  of  the  lift  fan  engine: 
the  powerplant,  the  transmission  shafts,  and  the  lift  fan.  The  thrust  vectoring  nozzle  on 
the  powerplant  and  the  mechanical  gearing  that  connected  the  transmission  shaft  to  the 
lift  fan  were  not  addressed.  However,  after  an  examination  of  current  technology  the 
following  assumptions  were  made  about  the  unexamined  components.  First,  a  swiveling 
nozzle  would  direct  the  thrust  of  the  powerplant  vertically  without  loss  when  required  for 
hover.  Second,  the  rotational  speed  of  the  lift  fan  should  be  held  to  within  25%  of  the 
rotational  speed  of  the  low  pressure  (LP)  engine  spool  in  order  to  limit  the  size  of  the 
mechanical  coupling. 


C.  BASELINE  AIRCRAFT 

The  baseline  aircraft  was  designed  to  have  the  characteristics  shown  in  Table  1. 
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Table  1.  Aircraft  Characteristics 


Aircraft  Parameter 

Value 

Length  (ft) 

114 

Width  (ft) 

13.2 

Wing  Span  (ft) 

96.6 

Wing  Chord  at  Root  (ft) 

25.6 

Aircraft  Weight  with  Two  Lift  Fan  Engines  (lbs) 

82500 

Maximum  Cargo  Weight  (lbs) 

37500 

Maximum  Aircraft  Weight  -  With  Two  Engines,  Without  Fuel  (lbs) 

120000 

D,  PERFORMANCE  ESTIMATES 

1,  Background 

The  drag  and  fuel  figures  for  the  baseline  aircraft  were  calculated  through  the  use 
of  the  AEDsys  software  program  (Mattingly  and  Pratt).  The  software  program  required 
initial  estimates  for  the  weight  of  the  aircraft  (including  fuel)  and  the  thrust  of  the  aircraft 
at  sea  level.  For  the  software  calculations,  the  total  aircraft  weight  (including  fuel)  was 
estimated  to  be  180000  lbs.  This  number  was  determined  by  adding  the  maximum 
aircraft  weight  of  120000  lbs  with  an  estimated  aircraft  fuel  weight  of  60000  lbs.  In 
addition,  the  thrust  was  estimated  to  be  60000  Ibf.  The  estimates  for  aircraft  weight  and 
thrust  were  then  used  to  calculate  the  drag  and  fuel  used  by  the  aircraft  over  the  following 
flight  profile:  Vertical  takeoff,  conventional  climb  to  20000  feet,  cruise  leg  at  20000  feet 
and  M=0.6,  and  then,  loiter  for  15  minutes  before  landing  vertically.  The  total  distance 
covered  during  the  climb  and  cruise  was  300  nautical  miles. 

2,  Results 

The  data  for  the  flight  profile  are  shown  in  Table  2.  The  data  showed  that  28000 
lbs  of  fuel  were  required  to  complete  one  300  nautical  mile  leg.  The  mission  fuel 
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requirement  (600  nautieal  miles  total)  was  found  on  the  assumption  that  the  mission 
would  involve  two  identieal  300  nautieal  mile  legs  with  a  landing  after  eaeh  leg.  Thus, 
the  total  mission  fuel  used  was  estimated  to  be  56,000  lbs.  It  was  noted  that  the  aireraft 
had  a  diffieult  time  aeeelerating  as  shown  by  the  large  amount  of  fuel  used  in  the  elimb. 
The  fuel  numbers  indieated  that  although  the  aireraft  funetioned  with  a  thrust  of  60,000 
lbs,  a  greater  amount  of  thrust  was  desired. 


Table  2.  Aireraft  Thrust  and  Fuel  Requirements  (Aireraft  Weight:  180000  lbs) 


Flight  Profile 

Drag  (lbs) 

Fuel  (lbs) 

Start  and  Hover  for  2.5  min 

175000 

3000 

Climb  to  20000  feet  (1) 

-22500 

11000 

Cruise  at  20000  feet  /  M=0.6  (1) 

19500 

7500 

Loiter  for  15  min  at  5000  feet 

15000 

3500 

Hover  for  2.5  min  and  Land 

175000 

3000 

Total  Fuel  for  One  300  Nautical  Mile  Leg 

28000 

Note  (1):  The  total  distance  covered  during  the  climb  and  cruise  is  300  nautical  miles 

E.  ENGINE  PERFORMANCE  CALCULATIONS 

The  program  used  to  ealeulate  the  engine  performanee  was  GASTURB  Version  9 
(Copyright  ©  2001,  Joachim  Kuzke).  The  ealeulations  were  made  using  the  Two  Spool  Mixed 
Flow  Turbofan  engine  eonfiguration  in  GASTURB.  Figure  1  is  a  generie  diagram  of  a 
Two  Spool  Mixed  Flow  Turbofan  engine  eonfiguration  and  it  gives  the  notation  for 
stations  within  the  engine. 
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d  Handling  bleed 


MFanLBCD.WMF  GasTurb 

Figure  1.  Two  Spool  Mixed  Flow  Turbofan  Engine  Configuration 

F.  INSTALLED  ENGINE  PERFORMANCE 

The  software  programs  used  in  the  present  projeet  calculated  the  uninstalled  thrust 
from  the  lift  fan  and  the  engine.  The  programs  did  not  take  into  account  the  loss  of  thrust 
once  the  engine  and  lift  fans  were  installed  into  the  aircraft.  Since  there  was  no 
established  way  to  calculate  the  installed  thrust,  it  was  assumed  that  there  was  a  10% 
difference  between  the  uninstalled  thrust  and  the  installed  thrust.  For  all  calculations  of 
thrust,  the  uninstalled  thrust  was  reduced  by  10%  to  provide  a  value  for  the  installed 
thrust  of  the  aircraft. 
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II.  POWERPLANT 


A.  DESIGN  CONSIDERATIONS 

There  were  two  options  when  determining  the  powerplant  for  a  lift  fan  engine: 
alter  an  existing  engine  or  design  a  new  one.  The  present  study  foeused  on  altering  an 
existing  engine  beeause  it  was  viewed  as  the  most  realistie  approaeh.  Although  the 
engine  would  not  provide  the  optimum  design,  it  would  result  in  a  system  more  likely  to 
reaeh  serviee  in  the  near  future. 

The  engine  for  the  F-35C  is  currently  the  only  lift  fan  that  is  approaching 
operational  status.  The  engine  is  a  low  bypass  turbofan  engine  that  has  provisions  for 
driving  the  lift  fan  from  the  Low  Pressure  (LP)  engine  spool.  Equation  1  is  the 
expression  of  the  balance  of  power  between  the  turbine  and  the  core  compressor,  bypass 
fan  and  lift  fan. 


(1  +  a) 


a 


(1) 


Mo 


where  f=fuel  to  air  ratio,  a=bypass  ratio,  Cp=ratio  of  specific  heats  in  the  turbine  (t)  and 
compressor  (c),  Tt=stagnation  temperatures  at  the  stations  depicted  in  Figure  1,  Pi^Power 

to  the  lift  fan,  and  mo  =engine  air  flow  rate.  If  the  turbine  is  designed  for  the  engine 
without  a  lift  fan,  it  must  operate  far  from  its  design  point  when  the  lift  fan  is  engaged.  If 
it  is  designed  for  operation  with  the  lift  fan  engaged,  it  will  be  far  off  design  when  the  lift 
fan  is  disengaged. 

In  conducting  the  present  study,  the  low  bypass  ratio  (a)  of  the  current  engine  was 
accepted,  but  a  larger  fraction  of  the  turbine  power  was  required  for  the  lift  fan.  This 
required  proposing  a  modification  of  the  turbine,  and  changing  the  nominal  design  point 
to  be  in  hover. 
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B,  BASELINE  ENGINE 

The  baseline  engine  used  to  begin  the  study  was  one  ereated  to  approximate  the 
engine  in  the  F-35C.  The  engine  specified  had  a  low  bypass  ratio  and  high  performance 
characteristics.  The  F-35C  engine  specifications  were  not  available,  but  representative 
specifications  for  the  technology  level  of  the  engine  were  estimated  as  shown  in  Table  3. 
The  complete  list  of  engine  parameters  entered  into  the  GASTURB  program  is  given  in 
Appendix  A. 


Table  3.  Estimated  Base 

ine  Engine  Parameters 

Engine  Parameter 

Value 

Design  Point 

30000  feet  and  M=L47 

Bypass  Ratio 

0.30 

Inlet  Corrected  Mass  Flow  (Ibm/sec) 

370 

Inner  Fan  Pressure  Ratio 

5 

Outer  Fan  Pressure  Ratio 

5 

High  Pressure  Compressor  Pressure  Ratio 

6.4 

Maximum  Burner  Temperature  (°R) 

3760 

Low  Pressure  Turbine  Pressure  Ratio 

2.017 

Ag  (in^) 

444.39 

Thrust  at  Design  Point  (Ibf) 

23599 

Length 

20  feet 

Hub  to  tip  ratio 

0.40 

C.  CHANGES  TO  THE  BASELINE  ENGINE 

The  baseline  engine  had  a  pressure  ratio  of  1.8  across  the  LP  turbine.  Since  the 
pressure  ratio  needed  to  be  increased  in  order  to  increase  the  work  output  of  the  turbine,  a 
redesign  of  the  turbine  was  proposed.  Other  than  the  modification  of  the  turbine,  all  the 
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other  turbo-machinery  components  of  the  baseline  engine  were  incorporated  into  the  new 
engine. 

Also,  the  baseline  engine  was  intended  for  a  high  performance  aircraft  with  a 
nominal  design  point  at  30000  feet  and  M=1.47.  Because  the  main  focus  of  the  heavy 
transport  aircraft  was  a  vertical  takeoff,  the  derivative  engine  was  designed  to  achieve 
maximum  performance  at  0  feet  and  M=0.0.  In  addition,  the  engine  was  optimized  for 
hover  operations  where  work  was  being  drawn  off  the  turbine  to  power  the  lift  fan. 
Although  this  meant  that  the  engine  would  run  inefficiently  in  cruise,  the  tradeoff  was 
needed  to  maximize  the  hover  performance. 


D,  DERIVATIVE  ENGINE 

1.  Referencing  Engine  Performance  to  New  Design  Point 

The  first  step  in  the  redesign  process  was  to  re-reference  the  performance  of  the 
engine  from  its  original  design  point  to  its  new  design  point  using  the  software  program 
GASTURB.  The  baseline  engine,  with  a  design  point  at  30,000  ft  and  M=1.47,  was  run 
off  design  at  0  feet  and  M=0.0.  During  the  run  the  corrected  high  pressure  spool  speed 
and  the  corrected  low  pressure  spool  speed  were  limited  to  100%.  The  off-design 
performance  parameters  at  0  feet  and  M=0.0  are  listed  in  the  first  column  in  Table  4. 
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Table  4.  Estimated  Lift  Engine  Parameters 


Engine  Parameter 

Baseline  Engine 

Derivative  Engine  at 

Off  Design 

Design  Point 

Flight  Condition 

0  Feet,  M=0.0 

0  Feet,  M=0.0 

Bypass  Ratio 

0.308 

0.308 

Inlet  Corrected  Mass  Flow  (Ibm/sec) 

369.635 

369.635 

Inner  Fan  Pressure  Ratio 

4.966 

4.966 

Outer  Fan  Pressure  Ratio 

4.966 

4.966 

Fligh  Pressure  Compressor  Pressure 

6.4 

6.4 

Ratio 

Maximum  Burner  Temperature  (°F) 

3347 

3760 

Low  Pressure  Turbine  Pressure  Ratio 

2.016 

1.818 

As  (in^) 

444.39 

454 

Thrust  (Ibf) 

33767 

37402 

LP  Shaft  Florse-power  (SFIP) 

42533 

42546 

The  basie  parameters  governing  the  inner  fan,  outer  fan  and  high  pressure  eompressor 
were  taken  from  the  off-design  ease  and  used  in  the  new  design  point  as  shown  in  Table 
4.  These  numbers  ensured  that  the  all  of  the  eomponents  exeept  the  turbine  operated  in 
the  same  manner.  This  proeess  left  the  turbine  as  the  only  turbo-maehinery  eomponent  to 
be  redesigned. 

2,  Turbine  and  Nozzle  Changes 

The  seeond  step  in  the  redesign  was  establishing  a  proeess  that  altered  the  engine 
to  produee  exeess  work  out  of  the  turbine  to  power  the  lift  fan.  The  eomponent  of  the 
engine  most  responsible  for  the  eontrol  of  the  turbine  work  output  was  the  area  at  station 
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8  (Ag)  (see  Figure  1).  Inereasing  Ag  inereased  the  pressure  drop  aeross  the  turbine  and 
resulted  in  an  inerease  in  work  produeed  by  the  turbine.  As  long  as  the  inereased  work 
was  drawn  off  to  power  the  lift  fan,  all  other  eomponents  of  the  engine  functioned  as 
before.  The  process  of  altering  Ag  and  controlling  the  turbine  work  output  was 
accomplished  using  the  LP  spool  mechanical  efficiency  parameter  in  the  GASTURB 
program.  Although  the  LP  spool  mechanical  efficiency  parameter  did  not  affect  a 
physical  engine  component,  reducing  the  parameter  properly  simulated  the  extraction  of 
shaft  power  for  a  lift  fan. 


3,  Engine  Changes  with  Changes  in  LP  Mechanical  Efficiency 

Table  5  shows  how  the  parameters  of  the  engine  varied  with  a  change  in  the  LP 
spool  mechanical  efficiency  (e).  Figures  2  and  3  show  the  variation  of  engine  thrust 
(FN),  and  total  SHP  with  the  size  of  Ag.  As  Ag  was  increased  the  Low  Pressure  Turbine 
Pressure  Ratio  increased.  The  increase  in  the  pressure  ratio  resulted  in  an  increase  in  the 
total  SHP  produced  by  the  LP  Spool.  Since  only  44440  SHP  was  required  to  drive  the 
engine  fan,  the  increased  work  was  removed  to  provide  power  to  the  lift  fan,  as  shown  by 
the  available  SHP.  The  figures  show  that  an  increase  in  A8  resulted  in  a  small  loss  of 
engine  thrust,  but  a  significant  increase  in  the  SHP  available  for  the  lift  fan. 
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Table  5.  Variation  of  Engine  Parameters  with  LP  Mechanical  Efficiency 


LP 

Mechanical 

Efficiency  (e) 

Ag  (in^) 

Thrust 

(Ibf) 

Low  Pressure 

Turbine  Pressure 

Ratio 

Total  LP 

SHP 

LP  SHP 

Available  for 

Lift  Lan 

1.00 

454 

37402 

1.818 

42546 

0 

0.95 

456 

37191 

1.880 

44786 

2239 

0.90 

459 

36952 

1.953 

47274 

4727 

0.85 

462 

36677 

2.038 

50054 

7508 

0.80 

467 

36358 

2.139 

53183 

10637 

0.75 

473 

35981 

2.261 

56728 

14182 

0.70 

482 

35524 

2.411 

60780 

18234 

0.65 

499 

34929 

2.600 

65435 

22910 

0.60 

536 

34099 

2.844 

70910 

28364 

0.55 

587 

33022 

3.169 

77357 

34811 

0.50 

658 

31598 

3.622 

85092 

42546 
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X  10 


Plot  of  AS  vs  Thrust  (FN) 


S  3.4 


AS  (in^) 


Figure  2.  Variation  of  Engine  Thrust  with  A8 


Figure  3.  Variation  of  Total  LP  SHP  with  A8 
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4,  Selection  of  the  Derivative  Engine  Design  Point 

Because  the  designed  called  for  the  thrust  produced  in  hover  to  be  fairly  well 
balanced  between  the  lift  fan  and  the  engine,  a  mechanical  efficiency  of  0.6  was  chosen 
for  the  design  point  of  the  engine.  The  final  design  parameters  of  the  engine  are  those  in 
Table  6. 


Table  6.  Estimated  Li 

f  Engine  Parameters 

Engine  Parameter 

Design  Point 

0  feet  and  M=0.0 

Bypass  Ratio 

0.308 

Inlet  Correeted  Mass  Flow  (Ibm/see) 

369.635 

Inner  Fan  Pressure  Ratio 

4.966 

Outer  Fan  Pressure  Ratio 

4.966 

High  Pressure  Compressor  Pressure  Ratio 

6.4 

Maximum  Burner  Temperature  (°F) 

3760 

Low  Pressure  Turbine  Pressure  Ratio 

2.844 

As  (in^) 

536 

Thrust  at  Design  Point  (Ibf) 

34099 

LP  Shaft  Horse-power  (SHP) 

70910 

LP  Shaft  Horse-power  (SHP)  Available  for  Lift  Fan 

28364 

Length 

20  feet 

Hub  to  tip  ratio 

0.40 
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E. 


ENGINE  PERFORMANCE  IN  CRUISE 


1,  Design  Considerations 

As  stated  earlier,  the  engine  was  optimized  for  the  hover  case  in  that  it  was 
designed  for  zero  altitude,  M=0.0  with  a  LP  spool  power  take-off  of  40%  for  the  lift  fan. 
For  all  other  operating  points,  the  engine  was  operated  off-design.  In  cruise,  because  the 
lift  fan  no  longer  drew  power  from  the  engine,  the  LP  spool  mechanical  efficiency  was 
required  to  be  set  at  99%.  The  increased  mechanical  efficiency  altered  the  engine 
performance  since  the  work  that  the  turbine  was  designed  to  produce  was  no  longer 
required.  The  additional  work  was  absorbed  by  the  bypass  fan.  This  required  that  the 
Low  Pressure  Spool  Speed  in  GASTURB  be  limited  to  ensure  that  the  engine  fan  and  the 
LP  Turbine  operated  at  acceptable  points  on  the  associated  component  maps,  and  always 
below  the  design  spool  speed.  In  addition,  the  nozzle  area  was  required  to  be  adjusted  at 
each  off-design  point  to  ensure  that  the  work  balance  between  the  engine  fan  and  the  LP 
turbine  was  such  that  engine  fan  provided  a  high  pressure  ratio  at  a  high  efficiency. 
Figure  4  and  5  show  the  input  screens  in  GASTURB  that  allowed  the  variation  of  the 
Low  Pressure  Spool  Speed  and  the  nozzle  to  be  controlled.  The  final  choices  of  Low 
Pressure  Spool  Speed  and  nozzle  area  were  based  on  selecting  the  combination  that 
provided  the  greatest  thrust  from  the  engine. 


Figure  4.  GASTURB  Limiters  Selection  Screen 
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Figure  5.  GASTURB  Modifiers  Seleetion  Sereen 


2,  Cruise  Performance 

Table  7  shows  the  performance  of  one  derivative  engine  at  three  forward  thrust 
operating  conditions.  Figure  6  and  7  show  component  maps  for  the  engine  fan  and  LP 
turbine  respectively  with  the  operating  condition  at  20000  ft  and  M=0.6.  The  remainder 
of  the  engine  component  maps  and  operating  conditions  for  the  points  in  Table  7  are 
given  in  Appendix  B.  The  data  indicated  that  the  derivative  engine  had  enough  thrust 
such  that  only  two  engines  were  required  to  operate  the  aircraft  up  to  approximately 
20000  feet.  The  small  difference  between  the  thrust  available  and  the  thrust  required, 
though,  indicated  that  aircraft  would  have  a  difficult  time  accelerating.  If  the  aircraft  was 
required  to  have  greater  performance  or  be  operated  above  20000  feet  then  the  aircraft 
would  require  additional  thrust.  It  can  be  seen  in  Figure  6  that  the  LP  compressor 
operating  point  was  relatively  near  the  stall  line  of  the  compressor.  This  was  considered 
acceptable  because  the  mission  of  the  aircraft  would  not  require  the  engine  to  operate 
under  flight  conditions  where  the  throttle  would  be  changed  rapidly. 
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Table  7.  Cruise  Performance  for  One  Engine 


Cruise 

Condition 

Change 
in  Nozzle 
Area  (%) 

Maximum 
LP  Spool 
Speed 
(rpm) 

Unistalled 

Thrust 

Available 

(Ibf) 

Installed 

Thrust 

Available 

(Ibf) 

Total 

Installed 

Thrust 

Available 

(Ibf) 

Thrust 

Required 

for 

Aircraft 

(Ibf) 

Oft 

M=0.0 

-2 

86.9 

20540 

18490 

36980 

N/A 

5000  ft 
M=0.55 

-14 

92.7 

17972 

16172 

32344 

25400 

20000  ft 
M=0.6 

-18 

91.8 

12300 

11070 

22140 

19500 

15 


Figure  7.  Engine  LP  Turbine  Map 


F.  BLEED  AIR 

The  baseling  engine  had  an  Ag  of  444  in  ,  as  shown  in  Table  3.  The  derivative  lift 
fan  engine  had  an  Ag  of  536  in^  as  shown  in  Table  6.  The  differenee  of  92  in^  represented 
a  signifieant  inerease  over  the  original  design.  Although  it  was  not  determined  to  be 
neeessary  for  the  present  design,  there  was  an  alternate  way  to  aehieve  the  desired  work 
out  of  the  engine  without  using  as  large  an  inerease  in  Ag.  The  method  involved  taking 
the  bypass  air  from  the  fan  and  bleeding  it  overboard.  In  GASTURB,  the  parameter  was 
ealled  the  Relative  Fan  Overboard  Bleed  (W_Bld/W13).  Table  8  showed  an  alternative 
derivative  engine  design  that  utilized  bypass  air  being  bled  overboard.  The  bleed  air 
method  worked  beeause  it  redueed  the  mass  flow  at  Ag  for  the  same  turbine  pressure  ratio 
and  flow  rate.  Thus,  Ag  was  not  required  to  be  as  large.  The  information  in  Table  8 
showed  that  if  99%  of  the  bypass  air  was  bled  overboard  through  an  auxiliary  nozzle,  the 
required  Ag  was  only  457  inehes  with  the  same  LP  shaft-horsepower  generated.  In 
addition,  the  air  bled  overboard  eould  be  used  as  additional  thrust,  or  for  direetional 
eontrol.  Assuming  that  the  auxiliary  nozzle  was  eonvergent.  Table  8  showed  that  the 
overboard  bypass  air  generated  2717  Ibf  of  thrust. 
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Table  8.  Estimated  Derivative  Engine  Parameters 


Engine  Parameter 

Design  Point  -  Single 

Thrust  Nozzle 

Design  Point  with  Fan 

Bleed  Air  Nozzle 

Design  Point 

0  feet  and  M=0.0 

0  Feet,  M=0.0 

Bypass  Ratio 

0.308 

0.308 

Inlet  Corrected  Mass  Flow  (Ibm/sec) 

369.635 

369.635 

Inner  Fan  Pressure  Ratio 

4.966 

4.966 

Outer  Fan  Pressure  Ratio 

4.966 

4.966 

High  Pressure  Compressor  Pressure 

Ratio 

6.4 

6.4 

Maximum  Burner  Temperature  (°F) 

3760 

3760 

Low  Pressure  Turbine  Pressure  Ratio 

2.844 

1.818 

Relative  Fan  Overboard  Bleed  (W13) 

0% 

99% 

As  (in^) 

536 

457 

Thrust  at  Design  Point  (Ibf) 

34099 

29161 

LP  Shaft  Horse-power  (SHP) 

70910 

70910 

Mass  Flow  from  Overboard  Bleed 

(Ibm/s) 

0 

78.3 

Thrust  from  Overboard  Bleed  (Ibf) 

0 

2717 

Total  Unistalled  Thrust  from  Engine  (Ibf) 

34099 

31878 

Total  Installed  Thrust  from  Engine  (Ibf) 

30700 

28700 
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III.  TRANSMISSION  SHAFTS 


A.  DESIGN  CONSIDERATIONS 

A  transmission  shaft  connected  the  powerplant  to  the  lift  fan.  The  design  of  the 
transmission  shaft  was  important  because  lift  fan  engines  in  general  required  transmission 
of  relatively  large  amounts  of  power  from  the  powerplant.  In  addition,  the  transmission 
shaft’s  length,  thickness  and  weight  had  an  impact  on  the  overall  design  of  the  aircraft. 
Thus,  the  transmission  shaft  design  was  an  integral  part  when  designing  a  lift  fan  engine. 


B,  FAN  AND  SHAFT  ROTATIONAL  SPEED 

Because  the  transmission  shaft  was  an  extension  of  the  Low  Pressure  (LP)  engine 
spool,  the  design  of  the  shaft  was  dictated  by  the  rotational  speed  of  the  LP  engine  spool. 
The  determining  factor  in  the  LP  engine  spool  speed  was  the  maximum  centrifugal 
tensile  stress  of  the  fan  blade.  The  stress  was  based  primarily  on  the  material  of  the 
blade,  the  flowpath  throughflow  area  and  the  blade’s  centrifugal  acceleration  (Mattingly, 

Heiser  and  Pratt,  2002).  The  relationship  between  these  factors  is  shown  in  equation  2 


AN^ 


3600O-, 


(2) 


2 

where  A  is  the  flowpath  throughflow  area  in  in  ,  N  is  the  rotational  speed  in  rpm,  Qc  is 
the  maximum  allowable  tensile  stress,  p  is  the  density  of  the  material  and  At/ Ah  is  the 
taper  ratio  of  the  blade  (assumed  to  be  0.8).  Based  on  current  titanium  blade  technology 
the  value  for  Gc  was  estimated  to  be  70000  ksi  and  the  value  for  p  was  taken  to  be  9.08 

slug/ft  .  Using  these  values  in  equation  2  resulted  in  a  value  of  1.00  x  10  in  -rpm  for 

2 

AN  .  The  throughflow  area  was  calculated  using 

A  =  (3) 

/ 


where  Rt  was  the  radius  of  the  fan  blade  at  the  tip  and  Rh/Rt  was  the  ratio  of  the  hub 
radius  to  tip  radius.  Examining  current  technology,  representative  values  for  Rt  and  Rh/Rt 
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were  39  inehes  and  0.4  respeetively.  Inserting  these  values  into  equations  2  and  3 
resulted  in  a  ealculation  of  10,000  rpm  for  the  rotation  speed  of  the  engine. 


C.  TRANSMISSION  SHAFT  THICKNESS 


1.  Governing  Equations 

The  rotational  speed  of  the  LP  spool  and  the  shaft  horsepower  transmitted  to  the 
lift  fan  were  used  to  determine  the  diameter  of  the  shaft  and  the  thiekness  of  the  shaft 
wall.  The  rotational  speed  (N)  and  the  shaft  horsepower  (SHP)  were  related  to  the  torque 
(T)  generated  by  the  shaft  by 


SHP*  63025 
N 


(4) 


The  torque  was  then  related  to  the  shear  stress  (t)  generated  on  the  shaft  with  a  thin 
walled  elosed  eireular  seetion  using  (Allen  and  Haisler,  1985) 


r  = 


(5) 


where  ro  is  the  outer  radius  of  the  eireular  seetion  and  ri  is  the  inner  radius  of  the  eireular 
eross  seetion.  Combining  the  two  equations,  the  shear  stress  in  the  shaft  was  related  to 
the  dimensions  of  the  shaft,  the  shaft-horsepower  transmitted  by  the  shaft  and  the 
rotational  speed  of  the  shaft,  as 


\26050*  SHP* 
Nn{r^^-r^) 


(6) 


Sinee  the  shaft-horsepower  and  the  rotational  speed  were  determined  in  the  design  of  the 
powerplant,  the  size  of  the  shaft  (and  thiekness)  were  required  to  be  large  enough  to 
ensure  that  the  maximum  allowable  shear  stress  in  the  shaft  was  not  exeeeded. 
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2,  Shaft  Design 

Since  the  transmission  shaft  was  an  extension  of  the  LP  engine  spool,  the  shaft 
was  modeled  with  similar  dimensions  and  materials  as  the  LP  engine  spool.  Based  on 
eurrent  teehnology  for  an  engine  with  a  39  ineh  fan  diameter  and  a  fan  blade  hub-to-tip 
ratio  of  0.4,  the  derivative  engine  was  estimated  to  have  a  LP  drive  shaft  diameter  of  5 
inehes  and  a  shaft  wall  thickness  of  0.25  inehes.  A  titanium  alloy  (TI-6A1-4V)  was 
assumed.  The  properties  of  the  alloy  are  shown  in  Table  9. 


Table  9.  Properties  of  Titanium  Alloy  (Ti-6A1-4V) 


Material 

Modulus  of  Elastieity 

(E)  (Mpsi) 

Unit  Weight  (w) 
(Ib/in^) 

Yield  Shear  Stress  (r) 

(ksi) 

Titanium 

Alloy 

(Ti-6A1-4V) 

17.4 

.16 

87 

Using  equation  6  with  the  values  shown  in  Table  10,  the  maximum  shear  stress  for  the 
transmission  shaft  was  ealeulated  to  be  21  ksi. 


Table  10.  Values  Used  for  Caleulation  of  Maximum  Shear  Stress 


Outer  Radius  (ro) 

Inner  Radius  (q) 

SHP 

N  (rpm) 

Maximum  Shear  Stress  (r) 

(in) 

(in) 

(ksi) 

2.5 

2.25 

28300 

10000 

21 

The  maximum  shear  stress  of  21  ksi  was  well  below  the  maximum  yield  shear  stress. 
The  ealeulations  indieated  that  the  size  of  the  transmission  shaft  was  appropriate  for  the 
power  transmitted  to  the  lift  fan.  It  was  noted,  however,  that  the  baseline  engine  LP  drive 
shaft  was  designed  to  transmit  42,500  shp.  The  additional  28,300  shp  supplied  to  the 
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lifting  fan  was  a  40%  increase  from  the  original  design.  These  numbers  indieated  that  the 
thickness  of  the  derivative  engine  drive  shaft  would  have  to  be  inereased  slightly  to 
eompensate  for  the  inereased  power. 


C.  TRANSMISSION  SHAFT  LENGTH 

1.  Critical  Speed 

The  length  of  the  transmission  shaft  was  defined  as  the  length  from  the  fan  on  the 
powerplant  to  the  gearbox  loeated  on  the  lift  fan.  The  primary  consideration  in  designing 
the  length  of  the  transmission  shaft  was  the  critical  speed  of  the  shaft.  Because  the  shaft 
was  rotating,  eentrifugal  forees  aeted  on  the  shaft.  A  perturbation  of  the  centrifugal  force 
eould  be  resolved  into  horizontal  and  vertical  components  that  would  trigger  vertieal  and 
horizontal  vibrations  of  the  shaft.  The  vibrations  become  violent  when  the  angular  speed 
of  the  shaft  eoineides  with  the  natural  frequency  of  the  non-rotating  shaft  (Den  Hartog,  J.P., 
1985).  This  angular  speed  is  referred  to  as  the  critical  speed. 


2.  Governing  Equations 

The  natural  frequency  of  the  rotating  shaft  (cOn)  depends  on  the  material  of  the 
shaft  and  the  length  of  the  shaft  (1)  as 


o)„ 


(V) 


where  an  is  a  numerieal  eonstant  which  varies  with  the  boundary  eonditions  of  the  shaft, 
E  is  the  Modulus  of  Elasticity,  I  is  the  Moment  of  Inertia,  and  p,  is  the  mass  per  unit 
length.  I  and  p  are  given  by 


/  = 


(8) 


and 
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(9) 


w* n{r^  -r^) 

~gl 


where  w  is  speeifie  weight  of  the  material  and  g  is  the  gravitational  eonstant  (386.09 
in/s  ).  Combining  equations  7,  8  and  9  and  inserting  the  value  for  g  gives 


=  a„ 


mMA* E{r^^  -r^) 

2773 


w{r^  -r.  )l 


(10) 


Sinee  the  rotational  speed  of  the  shaft  was  known,  equation  10  was  used  to  determine  the 
length  of  the  shaft  that  resulted  in  a  natural  frequeney  that  avoided  the  rotational  speed  of 
the  shaft. 


3.  Actual  Design 

Sinee  the  shaft  was  eonstrained  at  the  fan  of  the  powerplant  and  the  gearbox  of 
the  lift  fan,  the  shaft  was  modeled  initially  as  a  single  ‘elamped-elamped’  beam.  Solving 
the  Bemoulli-Euler  equation  for  a  beam  elamped  at  both  ends,  the  values  for  an  are  those 
shown  in  Table  11. 


Table  11.  Values] 

or  an  for  Clamped-Clamped  Beam 

Node 

1 

2 

3 

4 

Value  for  an 

22.373 

61.670 

120.912 

199.855 

Inserting  the  values  for  an  along  with  the  material  properties  of  the  alloy  and  plotting  the 
natural  frequeney  versus  length  resulted  in  the  graph  shown  in  Figure  8. 
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X  ■jO'*  Shaft  Length  vs  Natural  Frequency 


Figure  8.  Plot  of  Shaft  Length  Versus  Natural  Frequency 


The  figure  demonstrated  that  the  rotational  speed  of  10000  rpm  for  the  transmission  shaft 
approached  the  natural  frequency  of  a  ‘clamped-clamped’  shaft  at  lengths  of 
approximately  90  inches.  In  order  for  the  shaft  to  operate  effectively  it  would  have  to  be 
designed  to  be  shorter  or  longer  that  90  inches.  The  behavior  of  the  beam  at  the  second 
natural  frequency  was  not  considered  because  90  inches  was  felt  to  be  the  largest  length 
for  a  possible  shaft  design.  It  should  be  noted  that  although  there  were  many  options  for 
modeling  the  shaft  such  as  two  ‘clamped- free’  beams  coupled  together,  these  options 
were  not  considered  due  to  the  ability  of  the  ‘clamped-clamped’  beam  to  meet  the  desired 
design  goals. 
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IV.  LIFT  FAN 


A.  DESIGN  CONSIDERATIONS 

The  lift  fan  was  required  to  turn  the  shaft  horsepower  generated  by  the  powerplant 
into  useful  vertieal  thrust  for  the  aireraft.  Beeause  the  fan  was  mounted  horizontally,  it 
gave  rise  to  a  number  of  questions.  The  primary  one  was  whether  the  lift  fan  could 
produce  the  required  thrust.  The  second  question  concerned  the  dimensions  of  the 
resulting  fan  with  respect  to  the  limited  space  on  the  aircraft.  The  lift  fan  not  only  was 
required  to  provide  the  needed  thrust,  but  it  was  required  to  do  so  in  a  size  that  would  be 
compatible  with  the  aircraft. 


B,  AEROTHERMODYNAMIC  ANALYSIS 

The  first  step  in  the  design  was  an  aerothermodynamic  analysis  to  determine  the 
approximate  parameters  of  the  lift  fan.  The  analysis  used  one-dimensional  equations  to 
describe  the  performance  of  a  generic  fan.  In  order  to  calculate  lift  fan  thrust  (Tf),  the 
parameters  required  to  be  specified  were;  shaft  horsepower  supplied  to  the  lift  fan  (SHP), 
diameter  (D),  hub-to-tip  ratio  (Rh/Rq,  axial  mach  number  (Mai),  inlet  stagnation 
temperature  (Tn),  inlet  stagnation  pressure  (Pn),  inlet  pressure  loss  (iid),  nozzle  pressure 
loss  (Tin),  and  fan  polytropic  efficiency  (Cf).  The  sequence  of  equations  included  the 
expression  for  the  fan  face  area, 

4  =  (U) 

the  flow  rate  through  the  fan, 

PA  -1 

+  (12) 

^jyRT,,g,  2 

the  power  required  by  the  fan, 

778 

SHP  =  mfC^UT^-\)*—  (13) 
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the  relationship  between  fan  pressure  and  temperature  ratio, 


r 


/ 


(14) 


and  the  expression  for  the  speeific  thrust  if  the  fan  nozzle  is  correetly  expanded  to 
atmospheric  pressure. 


STf  = 


2r  RT,  TTf 


1 

/-I  gc 


yef  y-\ 

f 71, Tin)  "  -1] 


(15) 


Then,  the  total  thrust  from  the  lift  fan  is  given  by 


Tj.=mSTf  (16) 

Examining  current  technology  and  using  the  available  values  of  SHP  fixed  the 
parameters  shown  in  Table  12.  That  left  the  thrust  to  depend  on  the  parameters  of  D  and 
Mai.  It  was  assumed  that  there  were  two  engines,  with  each  lift  fan  required  to  have  a 
thrust  of  55000  Ibf  in  order  for  the  aircraft  to  hover.  Figure  9  shows  the  variation  of 
thrust  with  mass  flow.  Figure  10  shows  the  variation  of  mass  flow  with  D  and  Mai  at  an 
inlet  pressure  of  14.7  psi.  By  examining  the  two  figures  in  order,  it  was  seen  that  55000 
Ibf  required  a  mass  flow  of  4800  Ibm/sec.  A  mass  flow  of  4800  Ibm/sec  in  turn  required 
a  fan  diameter  of  12.7  feet  and  at  an  axial  Mach  number  (Mai)  of  0.6. 
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Table  12.  Initial  Lift  Fan  Parameters 


Engine  Parameter 

Design  Point 

0  feet  and  M=0.0 

Shaft  Horsepower 

28300 

Hub  to  Tip  Ratio  (Rh/Rt) 

0.3 

Inlet  Stagnation  Temperature  (°F) 

518.67 

Inlet  Stagnation  Pressure  (psi) 

14.696 

Inlet  Pressure  Ratio  (Ud) 

0.99 

Nozzle  Pressure  Ratio  (7i„) 

0.98 

Fan  Efficieney  (ef) 

0.9 
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Plot  of  Fan  Mass  Flow  vs  Fan  Diameter  for  28300  Shp 


C.  LIFT  FAN  PRELIMINARY  DESIGN 

1,  Design  Tool 

The  preceding  analysis  gave  the  overall  size  of  the  lift  fan.  A  more  in  depth 
analysis  was  required  to  determine  its  preliminary  design.  The  analysis  was  carried  out 
using  a  MATLAB  6.5  program  (Aaron,  Ryan,  2003).  A  listing  of  the  program  is  given  in 
Appendix  C.  The  program  required  the  following  inputs:  Fan  blade  tip  speed  (Ut),  hub- 

to-tip  ratio  (Rh/Rt),  mass  flow  (m),  inlet  stagnation  temperature  (Tn),  inlet  stagnation 
pressure  (Pti),  inlet  pressure  drop  (TTd),  nozzle  pressure  drop  (Tin),  inlet  relative  flow  angle 
on  the  meanline  (Pi),  rotor  solidity  on  the  meanline,  stator  solidity  on  the  meanline,  rotor 
diffusion  factor  on  the  meanline,  rotor  aspect  ratio  on  the  meanline  and  stator  aspect  ratio 
on  the  meanline.  From  these  inputs  the  program  calculated  the  size  of  the  fan,  the  thrust 
provided,  and  velocity  diagrams  at  hub,  mean  radius  and  tip.  In  using  the  code,  the  goal 
was  to  obtain  the  required  overall  thrust  while  not  exceeding  limits  for  the  diffusion 
factor,  at  any  radius,  in  either  the  rotor  or  stator  blade  row. 
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2. 


Determination  of  Initial  Lift  Fan  Paramaters 


a.  Fan  Blade  Tip  Speed  (Ud 

The  fan  blade  tip  speed  was  determined  from  the  allowable  centrifugal 

tensile  stress  of  the  fan  blade.  As  discussed  in  chapter  2,  paragraph  B,  the  stress 

depended  on  the  material  of  the  blade,  the  flowpath  throughflow  area  and  the  rotational 

speed.  Using  equation  4  and  assuming  the  blades  were  made  out  of  the  titanium,  a  value 

1.00  X  10  in  -rpm  was  calculated  for  AN  .  Using  equation  2  with  Rh/Rt  equal  to  0.3 

2 

and  converting  rpm  to  co  (rad/sec),  the  value  for  AN  was  written  as 

=266381 1.44  ft  V  (17) 

The  maximum  fan  blade  tip  speed  was  RtCO=1632  ft/s. 


b.  Inlet  Relative  Flow  Angle  ( Pi) 

The  inlet  relative  flow  angle  was  calculated  using  basic  blade  theory  as 
shown  in  Figure  1 1 . 
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The  value  of  P  at  the  tip  was  then 


_  ,  R.co  ^  ,  R,co 

y9  =  tan  '(^)  =  tan  '(- 


■) 


(18) 


With  T=518.67°F,  Mai=0.6,  and  RtCO=1632  ft/s,  the  maximum  P  was  ealeulated  to  be  68.2 
degrees. 


3.  Design  of  Lift  Fan 

a.  Design  Iteration 

The  design  of  the  Lift  Fan  involved  running  the  MATLAB  program 
several  times  while  altering  the  inlet  mass  flow  (mdot),  inlet  relative  flow  angle  on  the 
meanline  (Pi),  rotor  solidity  on  the  meanline,  stator  solidity  on  the  meanline,  rotor 
diffusion  faetor  on  the  meanline,  rotor  aspeet  ratio  on  the  meanline  and  stator  aspeet  ratio 
on  the  meanline.  The  purpose  of  eaeh  run  was  to  maximize  the  thrust  produeed  while 
eomp lying  with  the  eonstraints  listed  in  Table  13. 


Table  13.  Lift  Fan  Design  Constraints 


Engine  Parameter 

Design  Point 

0  feet  and  M=0.0 

Maximum  Shaft  Horsepower 

28300 

Maximum  Fan  Blade  Tip  Speed  (ft/s) 

1632 

Maximum  Inlet  Flow  Angle  on  the  Meanline  (Pi)  (deg) 

68.2 

Diffusion  Factor  at  Any  Point  on  Rotor  or  Stator 

0.6 

b.  Lift  Fan  Gearing 

Using  RtCO=1632  ft/s.  Table  14  showed  the  relationship  between  fan 
diameter  and  rotational  speed.  Previously,  the  LP  spool  of  the  engine  was  ealeulated  to 
have  a  rotational  speed  of  10000  rpm.  In  order  to  limit  the  size  of  the  reduetion  gear 
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between  the  engine  and  the  lift  fan,  the  lift  fan  was  desired  to  operate  at  a  minimum  of 
25%  of  the  LP  spool  rotational  speed  or  2500  rpm.  This  rotational  speed  equated  to  a 
desired  maximum  fan  diameter  of  12  feet. 


Table  14.  Diameter  vs  Rotational  Speed  (Fan  Blade  Tip  Speed=1632ft/s) 


Lift  Fan  Diameter  (ft) 

Rotational  Speed  (rpm) 

3 

10390 

6 

5195 

12 

2597 

18 

1732 

c.  Fan  Designs 

Two  potential  designs  for  the  lift  fan  are  shown  in  Table  15.  Detailed  fan 
information  is  given  in  Table  16.  The  first  design  had  a  diameter  of  12.1,  but  it  did  not 
produee  the  required  installed  thrust  of  55000  Ibf  The  seeond  design  produeed  the 
desired  amount  of  thrust  from  the  lift  fan,  but  it  had  a  diameter  of  18.4  ft  and  a  rotational 
speed  of  1247  rpm.  The  most  signifieant  parameter  in  the  ealeulations  was  the  fan  blade 
tip  speed.  A  high  fan  blade  tip  speed  was  required  for  smaller  fan  sizes  and  higher 
rotational  speeds.  The  drawbaek  to  the  high  fan  blade  tip  speed  was  a  reduetion  in  thrust. 
The  first  design  was  ehosen  as  the  final  design  beeause  it  met  the  desired  rotational 
speed.  In  addition,  the  first  design  had  less  depth  to  the  fan.  The  maximum  width  of 
both  the  rotor  and  the  stator  blades  was  6  inehes.  The  problem  with  the  seleetion  was  that 
additional  thrust  was  required  for  the  aireraft  to  hover. 
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Table  15.  Potential  Lift 

^an  Designs 

Engine  Parameter 

First  Design 

Seeond  Design 

Design  Point 

0  feet  and  M=0.0 

0  Feet,  M=0.0 

Hub  to  Tip  Ratio  (Rh/Rt) 

0.3 

0.3 

Inlet  Stagnation  Temperature  (°F) 

518.67 

518.67 

Inlet  Stagnation  Pressure  (psi) 

14.696 

14.696 

Inlet  Pressure  Ratio  (rid) 

0.99 

0.99 

Nozzle  Pressure  Ratio  (rin) 

0.98 

0.98 

Maximum  Fan  Blade  Tip  Speed  (ft/s) 

1632 

1200 

Inlet  Flow  Angle  on  the  Meanline  (Pi)  (deg) 

62 

62 

Diffusion  Faetor  on  the  Meanline 

.205 

.205 

Rotor  Solidity  on  the  Meanline 

0.42 

.44 

Stator  Solidity  on  the  Mealine 

0.43 

.44 

Rotor  Aspeet  Ratio  on  the  Meanline 

13 

13 

Stator  Aspeet  Ratio  on  the  Meanline 

13 

13 

Mass  Flow  (Ibm/s) 

3910 

7050 

Diameter  (ft) 

12.1 

18.4 

Rotational  Speed  (rpm) 

2576 

1247 

Shaft  Horsepower 

28005 

28041 

Uninstalled  Thrust  (Ibf) 

43581 

60977 

Installed  Thrust  (Ibf) 

39200 

54900 
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Table  16.  Detailed  Fan  Design  Parameters 


Engine  Parameter 

First  Design 

Seeond  Design 

Number  of  Rotor  Blades 

32 

34 

Rotor  Blade  Length  (in) 

50.8 

77.2 

Rotor  Blade  Width  at  Tip  (in) 

5.99 

8.97 

Number  of  Stator  Blades 

33 

34 

Stator  Blade  Size  at  Tip  (in) 

5.94 

8.97 

Rotor  Diffusion  Faetor 

At  Hub 

0.5967 

0.5921 

At  Meanline 

0.205 

0.205 

At  Tip 

0.0959 

0.1680 

Stator  Diffusion  Faetor 

At  Hub 

0.5762 

0.5800 

At  Meanline 

0.2627 

0.2645 

At  Tip 

0.1680 

0.1690 
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V.  FINAL  AIRCRAFT  DESIGN 


A.  ADDITION  OF  THRUST 

1.  Requirement 

The  previous  steps  left  the  derivative  lift  fan  engine  design  for  the  aireraft,  short 
of  its  target.  The  first  deficiency  was  the  thrust  required  in  hover.  A  summary  of  the 
thrust  calculated  in  the  previous  chapters  is  provided  in  Table  17. 


Table  17.  Summary  of  Thrust  Calculations 


Installed  Thrust  from 
Poweplants  (Ibf) 

Installed  Thrust  from 
Lift  Fans  (Ibf) 

Total  Thrust  (Ibf) 

Maximum  Weight 
of  Aircraft  (lbs) 

61400 

78400 

139800 

175000 

The  numbers  showed  that  an  additional  35200  Ibf  of  thrust  was  required  to  hover  the 
aircraft  at  maximum  weight.  The  second  deficiency  was  the  limited  capabilities  of  the 
engines  during  cruise.  The  small  margin  between  the  thrust  available  and  the  thrust 
required  at  20000  feet  and  M=0.6  (Table  7)  indicated  that  the  aircraft  would  have  a 
difficult  time  achieving  a  cruise  altitude  of  20000  feet. 


2,  Solution 

A  proposed  solution  for  the  two  deficiencies  of  the  derivative  engines  was  the 
addition  of  two  baseline  engines,  having  the  performance  given  in  Table  3.  A  summary 
of  the  additional  engine’s  performance  is  shown  in  Table  18. 


Table  18.  F-35C  Engine  Performance 


Operating  Point 

Uninstalled  Thrust  (Ibf) 

Installed  Thrust 

(Ibf) 

Total  Thrust  (Ibf)  for 
Two  Engines 

0  feet  and  M=0.0 

33767 

30400 

20000  feet  and  M=0.6 

21241 

19020 
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Assuming  that  each  engine  added  6000  lbs  of  weight  to  the  design,  the  addition  of  the 
baseline  engines  provided  the  thrust  required  to  hover  as  shown  in  Table  19. 


Table  19.  Summary  of  Thrust  Calculations 


Installed  Thrust  from 
Powerplants  (Ibf) 

Installed  Thrust  from 
Lift  Fans  (Ibf) 

Total  Thrust  (Ibf) 

Maximum  Weight 
of  Aircraft  (lbs) 

122200 

78400 

200600 

187000 

In  addition  to  the  extra  thrust  for  hover,  the  38040  Ibf  of  thrust  at  20000  feet  and  M=0.6 
was  significantly  greater  than  the  19500  lbs  of  thrust  required  (Table  7).  Thus,  the 
additional  engines  provided  an  improved  cruise  capability  at  high  altitudes  as  well  as  an 
increased  hover  capability. 


B.  DESIGN  OF  TRANSMISSION  SHAFT 

The  transmission  shaft  was  designed  for  the  shortest  length  between  the 
powerplant  and  the  lift  fan.  The  final  shaft  configuration  is  shown  in  Figure  12.  Segment 
A  was  the  part  of  the  shaft  that  connected  the  powerplant  and  the  reduction  gear.  It  was 
set  at  3  feet  to  provide  room  for  the  intake  for  the  powerplant  and  it  turned  at  10000  rpm. 
Segment  B  was  the  part  of  the  shaft  that  connected  the  reduction  gear  to  the  lift  fan.  It 
was  set  at  6.5  feet  to  provide  clearance  between  the  fan  and  the  reduction  gear.  In 
addition,  it  turned  at  2500  rpm.  Figure  8  showed  that  the  lengths  of  both  shafts  were 
satisfactory  for  avoiding  the  natural  frequency  of  the  shaft  when  they  were  turning  at 
their  designed  rotational  speed. 


Figure  12.  Lift  Fan  Engine  Installation  Dimensions 
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C.  FINAL  AIRCRAFT  DESIGN 

The  final  aircraft  design  was  powered  by  two  lift  fan  engines  and  two  baseline 
engines.  A  diagram  of  the  aircraft  is  shown  in  Figures  13  and  14.  The  important  things 
to  note  were  the  large  increase  in  the  wing  area  required  to  incorporate  the  lift  fan.  In 
addition,  the  engines  were  spaced  on  the  fuselage  in  order  to  balance  the  thrust  with 
respect  to  the  center  of  gravity  of  the  aircraft. 


Figure  13.  Top  View  of  Heavy  Lift  Aircraft 
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Figure  14.  Perspective  View  of  Heavy  Lift  Aircraft 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Current  technology  engines  can  be  used  in  the  design  of  lift  fan  engines  for  a 
heavy  lift  aircraft.  The  steps  for  designing  derivative  lift  fan  engines  with  otherwise 
current  technology  were  successfully  outlined  in  this  project.  The  steps  involved  creating 
a  derivative  engine  by  optimizing  a  baseline  engine  for  hover  conditions.  When  the 
derivative  engine  was  used  to  power  a  large  lift  fan,  the  resulting  lift  fan  engine  generated 
a  significant  amount  of  vertical  thrust.  Although  the  thrust  from  the  lift  fan  engine  was 
limited  by  aircraft  constraints,  it  was  still  significantly  more  thrust  than  could  be 
produced  from  a  baseline  engine.  The  addition  of  two  derivative  lift  fan  engines  to  two 
baseline  engines  provided  a  heavy  lift  aircraft  with  enough  thrust  to  take-off  vertically  as 
well  as  cruise  at  speeds  greater  than  M=0.6. 


B,  RECOMMENDATIONS 

1.  Perform  a  complete  mission  study  in  order  to  determine  precisely  the 
aircraft  weight  and  fuel  requirements. 

2.  Develop  a  tool  to  carry  out  the  preliminary  design  of  the  LP  turbine. 

3.  Develop  tools  that  provide  more  thorough  design  of  the  transmission  shaft 
and  reduction  gear. 

4.  Develop  more  tools  to  allow  the  designer  to  examine  the  aircraft  structural 
constraints  when  determining  the  size  of  the  lift  fan. 
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APPENDIX  A 


Mixed  Flow  Turbofan  Inputs 

Altitude  (ft)  30000 

Delta  T  from  ISA  (°R)  0 

Relative  Humidity  0 

Maeh  Number  1.47 

BASIC  DATA 

Intake  Pressure  Ratio  1 

Inner  Fan  Pressure  Ratio  5 

Booster  Map  Type  (0/1/2)  0 

Outer  Fan  Pressure  Ratio  5 

Compressor  Interduet  Pressure  Ratio  0.99 

HP  Compressor  Pressure  Ratio  6.4 

Bypass  Duet  Pressure  Ratio  0.97 

Turbine  Interduet  Pressure  Ratio  0.98 

Design  Bypass  Ratio  0.3 

Burner  Exit  Temperature  (°R)  3760 

Burner  Design  Effieieney  0.9995 

Burner  Partload  Constant  1 .6 

Fuel  Heating  Value  (BTU/lb)  1 8552.4 

Handling  Bleed  Loeation  0 

Overboard  Bleed  (Ib/s)  0 

Power  Offtake  (bp)  67.051 1 

HP  Spool  Meebanieal  Effieieney  1 

LP  Spool  Meebanieal  Effieieney  1 

Burner  Pressure  Ratio  0.97 


43 


Turbine  Exit  Duct  Pressure  Ratio  0.98 

Hot  Stream  Mixer  Pressure  Ratio  0.98 

Cold  Stream  Mixer  Pressure  Ratio  0.99 

Mixed  Stream  Pressure  Ratio  1 

Mixer  Efficiency  0.5 

Design  Mixer  Mach  Number  0.2 

Design  Mixer  Area  (in^)  0 

Nozzle  Thrust  Coefficient  1 

Design  Nozzle  Petal  Angle  (°)  25 

AIR  SYSTEM 

Relative  Handling  Bleed  to  Bypass  0 

Relative  Enthalpy  of  HP  Handling  Bleed  0 

Relative  HP  Leakage  to  Bypass  0 

Relative  Overboard  Bleed  (W_BLdAV25)  0.005 

Relative  Enthalpy  of  Overboard  Bleed  1 

NGV  Cooling  Air  (W_CL_NGVAV25)  0.05 

LPT  Cooling  Air  (W_CL/W25)  0.03 

Relative  Enthalpy  of  LPT  Cooling  Air  0.6 

HPT  Cooling  Air  (W_CLAV25)  0.05 

Relative  HP  Leakage  to  LPT  Exit  0 

Relative  Fan  Overboard  Bleed  (W_BldAV13)  0 

MASS  FLOW  INPUT 

Inlet  Corrected  Flow  (W2Rstd)  (Ib/s)  370 

LPC  Efficiency 

Isentropic  Inner  LPC  Efficiency  0.89 

Isentropic  Outer  LPC  Efficiency  0.88 

LPC  Design 
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Nominal  LP  Spool  Speed  10000 

HPC  Efficiency 

Isentropic  HPC  Efficiency  0.86 

HPC  Design 

Nominal  HP  Spool  Speed  18000 

HPT  Efficiency 

Isentropic  HPT  Efficiency  0.9 

LPT  Efficiency 

Isentropic  LPT  Efficiency  0.91 
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P13q2  LPC  Outer  Pressure  Ratio 


APPENDIX  B 


COMPONENT  MAPS  FOR  0  FEET  AND  M=0.0  (CRUISE) 

HPC 


LPT 
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1.05 


Pressure  Ratio  P45/P5  P13q2  LPC  Outer  Pressure  Ratio 


COMPONENT  MAPS  FOR  5000  FEET  AND  M=0.55  (CRUISE) 
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Pressure  Ratio  P45/P5  P13q2  LPC  Outer  Pressure  Ratio 


COMPONENT  MAPS  FOR  20000  FEET  AND  M=0.6  (CRUISE) 
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APPENDIX  C 


Lift  Fan  Code 


Developed  by  Prof.  Ray  Shreeve  and  converted  to  Matlab®  by  Ryan  Aaron 


%  IMPULSE  STAGE  COMPRESSOR  DESIGN  PROGRAM  "CDESIGN" 

%  UPDATED  1-27-03  (R. P . SHREEVE)  FOR  LIFT  FAN  THRUST  AND  POWER 

g, 

o 

% - DECLARE  INTEGER  VARIABLES 

% - III - INPUT  DATA  STATEMENTS - III - 

q, 

o 

% - III — NECESSARY  BOUNDARY  CONDITIONS  —  III - 

0, 

o 

Utip=I000;  %  WHEEL  SPEED  AT  THE  FIRST  ROTOR  TIP  (FT/SEC) 

Rhub_tip=. 3;  %  HUB-TO-TIP  RADIUS  RATIO  AT  THE  ROTOR  FACE 

Kblock=.98;  %  INLET  BLOCKAGE  FACTOR  (TO  ALLOW  FOR  HUB  AND 

g, 

o 

%SHROUD  BOUNDARY  LAYERS 


0, 

o 


Mdot=3250; 

Pt0=14.7; 

TtO=5I8.7; 


Pid=. 99; 
Pin=. 98; 

q, 

o 


Cp= .238; 
Gc=32 .174; 
J=778; 
Gamma=l . 4 ; 

q, 

o 

q, _ _ _ 

o 


%  FLOW  RATE  (LBS /SEC) 

%  STAG. PRESSURE  INTO  FIRST  ROTOR 
%  STAG. TEMP. INTO  FIRST  ROTOR 


///--Fan  inlet  and  nozzle  loss--/// 


III - NECESSARY  CONSTANTS - III - 

%  SPECIFIC  HEAT  AT  CONSTANT  PRESSURE 
%  GRAVITATIONAL  CONSTANT  (LBF/LBM. FT/SEC2 ) 
%  MECHANICAL  EQUIV.  OF  HEAT  (FT.LBS/BTU) 

%  RATIO  OF  SPECIFIC  HEATS 

III - SECONDARY  VARIABLES - III - 


g, 

o 

Phir=l; 

Phis=l ; 

Lr=0 ; 

Ls=0; 

Nr=0; 

Ns=0; 

q, 

o 

% - III  PRIMARY 

q, 

o 

Betlm=54 . 6*pi/180; 

q, 

o 

Sigrm=l .25; 

Sigsm=l .25; 


%  RATIO  OF  AXIAL  VEL.  OUT/ INTO  ROTOR 
%  RATIO  OF  AXIAL  VEL.  OUT/ INTO  STATOR 
%  COEFF.  OF  1/RBAR  TERM  IN  SOLIDITY  VARIATION 
%  FOR  ROTOR 
%  SAME  FOR  THE  STATOR 

%  COEFF.  OF  REAR  TERM  IN  SOLIDITY  VARIATION  FOR 
%  ROTOR 

%  SAME  FOR  THE  STATOR 

CHOICES  III - 

%  INLET  RELATIVE  FLOW  ANGLE  ON  THE  MEANLINE 

%  ROTOR  SOLIDITY  ON  THE  MEANLINE 
%  STATOR  SOLIDITY  ON  THE  MEANLINE 
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Drm= .4636; 
Arr=6; 

Ars=8 ; 
Clrhr= .006; 

Clrhs= .006; 

0, 

o 

o 

o 


% - INPUT  COMPLETE 

q, 

o 


%  ROTOR  DFACTOR  ON  THE  MEANLINE 
%  ROTOR  ASPECT  RATIO  ON  THE  MEANLINE 
%  STATOR  ASPECT  RATIO  ON  THE  MEANLINE 
%  ROTOR  TIP  CLEARANCE /BLADE  HEIGHT 
%  STATOR  TIP  CLEARANCE /BLADE  HEIGHT 


STAGE  CALCULATION  STARTS 


Ptl=PtO; 

TtI=TtO; 

q, 

o 

%  BEGIN  THE  CALCULATION  OF  EACH  NEW  STAGE  HERE 

q, 

o 

disp ( 'RESULTS  FOR  LIFT  FAN ' ) ; 

disp(' - '); 

disp ( '  ' ) ; 

disp('TIP  SPEED  =  ');  disp (Utip) ; disp (' Ft/Sec ') ; 
disp ( ' FLOW  RATE= ' ) ; disp (Mdot) ; disp ( ' Ibs/sec ' ) 
disp ( ' Pt= ' ) ; disp ( Ptl ) ; disp ( ' psia ' ) ; 
disp ( ' Tt= ' ) ; disp (Ttl ) ; disp ( ' Deg . R ' )  ; 
disp ( '  '  )  ; 

q, 

o 

% - START  CALCULATIONS 

q, 

o 

% - INLET  FLOW  CONDITIONS 

q, 

o 

Uin=Utip*  ( 1+Rhub  tip)/2; 

Vtl=sqrt (2*Cp*Gc*J*TtI) ; 

Xulm=Um/Vtl ; 

XIm=Xulm/tan (Betlm) ; 

MIm=FNMf romx (Xlm, Gamma) ; 

Mwlm=Mlm/cos (Betlm) ; 

Mulm=Mlm*tan (Betlm) ; 

Gammastarl=FNGammastar (Xlm, Gamma) ; 

Al=Mdot*Vtl/  (144*Ptl*Kblock;*Gc*Gammastarl)  ; 

Rm=sqrt (Al* (1+Rhub  tip) / (4*pi* (1-Rhub  tip) ) ) ; 

H=2*Rm* (1-Rhub  tip) / (1+Rhub  tip) ; 

Omega=Um/Rm; 

Phil=l/ tan (Betlm) ; 

disp ('MEAN  RADIUS  =  '); disp (Rm* 12 ); disp ( '  INCHES') 
disp ('BLADE  HEIGHT  =  '); disp (H* 12 ); disp ( '  INCHES'); 

q, 

o 

% - ACROSS  THE  ROTOR 

q, 

o 

Bet2m=FNInv  dfact (Drm, Phir, Sigrm, Betlm) ; 

Rstm=. 5* (1+ (Phir*Phil*tan (Bet2m) ) ) ; 
disp  (  '  '  )  ; 

disp ( 'VELOCITY  DIAGRAM:  FOR  BETA  1=  ');  disp (Betlm) ; disp  (  ' 

DFACTOR=') ; disp (Drm) ; 

disp ('BETA  2=  ' ) ; disp (Bet2m) ; 

Alphalm=0 ; 

Alpha2m=atan (2* (1-Rstm) / (Phir*Phil) ) ; 
disp ('ALPHA  2=  ' ) ; disp (Alpha2m) ; 
disp  (  '  '  )  ; 


ROTOR 
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disp ( ' REACTION=  ' ) ; disp (Rstm) ; 

X2m=Phil*Phir*Xulm/cos (Alpha2m) ; 

Tau=l+4* (1-Rstm) *Xulm^2; 

Y2m=X2m/ sqrt (Tau)  ; 

M2m=FNMf romx (Y2m, Gamma) ; 

Mw2m=M2m*cos (Alpha2m) /cos (Bet2m) ; 

g, 

o 

g, 

o 

g, 

o 

% - ACROSS  THE  STATOR 

0, 

o 

Y3m=Phis*Y2m; 

M3m=FNMf romx (Y3m, Gamma) ; 

Alpha3=0 ; 

Dsm=FNDfact (Phis, Sigsm, Alpha2m, Alpha3 ) ; 
disp ( ' STATOR  DFACTOR=' ) ;disp (Dsm) ; 
disp  (  '  '  ) 

g, 

o 

% - DECIDE  ON  BLADE  NUMBERS 

g, 

o 

Chordrm=H/Arr; 

Zr=2*pi*Rm*Sigrm/ Chordrm; 

disp ('NUMBER  OF  ROTOR  BLADES  IS  ' ) ; disp ( Zr ) ; 

Zr=input ( 'ENTER  WHOLE  NUMBER  REQUIRED  '  )  ; 
disp ('CHOSEN  WHOLE  NUMBER  IS  ' ) ; disp ( Zr ) ; 

Sprm=2  *pi*Rm/ Zr ; 

Chordrm=Sprm*Sigrm; 

Chordsm=H/Ars; 

Zs=2  *pi*Rm* Sigsm/ Chordsm; 

disp ('NUMBER  OF  STATOR  BLADES  IS  ' ) ; disp ( Zs ) ; 

Zs=input ( 'ENTER  WHOLE  NUMBER  REQUIRED  '); 
disp ('CHOSEN  WHOLE  NUMBER  IS  ' ) ; disp ( Zs ) ; 

Spsm=2  *pi*Rm/ Zs ; 

Chordsm=Spsm* Sigsm; 

Spsmh=Spsm/H; 

Sprmh=Sprm/H; 
disp  (  '  '  )  ; 

disp ( 'BLADING  SIZES:  FOR  ASPECT  RATIO  ROTOR= '); disp (Arr) ; disp 

STATOR=' ) ;disp (Ars) ; 
disp  (  '  '  ) 

disp ('ROTOR  CHORD  =  ');  disp (Chordrm* 12 ); disp  (  '  inches'); 
disp ('SPACE  =  '); disp ( Sprm* 12 ); disp ( '  '); 

disp ('STATOR  CHORD  =  '); disp (Chordsm* 12 ); disp ( '  inches'); 
disp ('SPACE  =  '); disp ( Spsm* 12 ); disp ( '  '); 

disp ( '  '  )  ; 

g, 

o 

% - SEC.  FLOW  AND  T.C.LOSS  FROM  MEAN  LINE  CONDITIONS 

g, 

o 

Omegarstc=FNOmegastc (Sigrm, Betlm, Bet2m, Sprmh,  Clrhr) ; 
Omegasstc=FNOmegastc (Sigsm, Alpha2m, Alpha3 , Spsmh, Clrhs ) ; 
disp (' SEC. &  T.C.  LOSS  FOR  ROTOR=  ');disp(Omegarstc); 
disp (' SEC. &  T.C.  LOSS  FOR  STATOR=  ' ) ; disp (Omegasstc) ; 
disp  (  '  '  )  ; 

disp(' - ') 

g, 

o 
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RADIAL  DISTRIBUTION  OF  CONDITIONS  &  PERFORMANCE 


0, 

o 

disp ( '  ' ) ; 

disp ( 'HUB-TO-TIP  CALCULATIONS ' ) ; 

disp(' - '); 

disp ( '  '  )  ; 

N=2;  %  N=NUMBER  OF  STREAMLINES 

disp ('FOR  SOLIDITY  VARIATION  IN  THE  ROTOR,  Lr= '); disp (Lr) ; disp ( ' 
Nr=' ) ;disp (Nr) ; 

disp ('AND  IN  THE  STATOR,  Ls= '); disp ( Ls ); disp ( '  and  Ns= ' ) ; disp (Ns ) ; 

disp(' - 

Invn=l /N; 

Rbart=2/ (1+Rhub  tip) ; 

Rbarh=Rbart*Rhub  tip; 

q, 

o 

% - SET  UP  COUNTER  TO  AVERAGE  P.R. 

q, 

o 

Number=0 ; 

Prttsum=0 ; 
for  I=0:Invn:I 
Number=Nuinber+l ; 

Rbar=l+ ( 1-Rhub  tip) * (2*1-1) / (1+Rhub  tip); 
disp ( ' AT  RBAR= ' ) ; di sp ( Rbar ) ; 

Xl=Xlm; 

Xul=Rbar*Xulm; 

Xwl=sqrt ( (X1^2+Xul^2)  )  ; 

Betl=asin (Xul /Xwl ) ; 

Ml=Mlm; 

Mwl=Ml /cos (Betl )  ; 

Alpha2=atan (tan (Alpha2m) /Rbar) ; 

Bet2=atan ( (Rbar/ (Phil*Phir) ) -tan (Alpha2m) /Rbar) ; 

Y2=Y2m*cos (Alpha2m) /cos (Alpha2) ; 

M2 =FNMfromx (Y2 , Gamma) ; 

Mz2=M2*cos (Alpha2)  ; 

Mw2=Mz2 / cos (Bet2 ) ; 

Yw2=FNXfromm (Mw2 , Gamma)  ; 

Y3=Y2*cos (Alpha2) *Phis; 

M3=FNMfromx (Y3 , Gamma) ; 

Sigr=Sigrm* (Lr+Rbar+Nr*Rbar^2) /Rbar;  %  VARIATION  IN  SOLIDITY 

Sigs=Sigsm* (Ls+Rbar+Ns*Rbar^2) /Rbar;  %  VARIATION  IN  SOLIDITY 

Dr=FNDf act ( Phir , Sigr , Betl ,  Bet2  )  ; 

Ds=FNDfact (Phis, Sigs , Alpha2 , Alpha3 ) ; 
disp ( '  '  )  ; 

disp ( 'VELOCITY  DIAGRAM:'); 
disp ('BETA  1  =  '); disp (Betl ) ; 
disp ('BETA  2  =  ');disp(Bet2); 
disp('ALPHA2  =  ' ) ; disp (Alpha2 ) ; 
disp  (  '  '  ) 

disp ('ROTOR  DFACTOR  =  ');disp(Dr); 
disp ('STATOR  DFACTOR  =  ');disp(Ds); 

q, 

o 

%  CALCULATE  BLADE  SIZES 

0, 

o 

Spr=Sprm*Rbar ; 

Sps=Spsm*Rbar ; 


and 
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Chordr=Spr*Sigr ; 

Chords=Sps*Sigs; 
disp ( '  '  )  ; 

dispC  BLADE  SIZES:'); 

disp ('ROTOR  CHORD  =  '); disp (Chordr* 12 ); disp (' ins .') ; 
disp (' SPACING  =  '); disp ( Spr* 12 ); disp (' ins .') ; 
disp (' SOLIDITY  =  '); disp ( Sign) ; 

disp ('STATOR  CHORD  =  '); disp (Chords* 12 ); disp (' ins .') ; 
disp (' SPACING  =  '); disp (Sps* 12) ; disp (' ins .') ; 
disp (' SOLIDITY  =  ' ) ; disp ( Sigs ) ; 
disp ( '  '  )  ; 

q, 

o 

g, 

o 

% - PERFORMANCE 

q, 

o 

Gogl=Gamma/ (Gamma-1) ; 

Tlttl  =  l-Xl''2; 

Plptl= (1-X1^2) ^Gogl; 

Trlttl=Tlttl+Xwl''2; 

Prlptl=Trlttl''Gogl; 

o 

o 

% - TEST  FOR  SHOCK  AND  SHOCK  LOSS 

g, 

o 

if  Mwl>l 

Sptratio=FNShock;ptratio  (Mwl ,  Gamma)  ; 

Prlptl=Prlptl*Spt ratio; 

Spratio=FNShock;pratio  (Mwl ,  Gamma)  ; 

Plptl=Plptl* Sprat io; 

disp ('A  SHOCK  WAS  FOUND  AT  MWl=');disp(Mwl); 

disp ( ' PR (T-T) =  ' ) ; disp ( Spt ratio) ; disp ( ' 

' ) ;disp (Spratio) ; 

end  %  AFTER  SHOCK  CONDITIONS  FROM  HERE 

Omegarp=FNOmegap (Dr, Sigr, Betl ,  Bet2 ) ; 
Omegart=Omegarp+Omegarstc; 

Pr2ptl=Prlptl-0megart* (Prlptl-Plptl) ; 

Tr2ttl=Trlttl; 

T2ttl=Tr2ttl-Yw2’^2; 

P2ptl= ( (T2ttl/Tr2ttl) ^Gogl) *Pr2ptl; 

Pt2ptl=P2ptl /  (  (T2ttl/Tau)  ^'Gogl)  ; 

q, 

o 

% - ACROSS  THE  STATOR 

q, 

o 

Omegasp=FNOmegap (Ds, Sigs , Alpha2 , AlphaS ) ; 

Omegas t=Omegasp+Omegasstc; 

Pt3ptl=Pt2ptl-Omegast* (Pt2ptl-P2ptl) ; 

Tt3ttl=Tau; 

disp ( ' PT3/PT1=  ' ) ;disp (Pt3ptl) ; 
disp ( '  '  )  ; 

disp(' - 

Prttsum=Pt3ptl+Prttsum; 

end 

Prttavg=Prttsum/Number ; 

disp ( 'AVERAGE  PRESSURE  RATIO (T-T) =  ' ) ; disp ( Prttavg) ; 

g, 

o 

% - CALCULATE  CONDITIONS  FOR  NEXT  STAGE 


PR(S-S) 
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Kblock3=. 98;  %BLOCKAGE  FACTOR 

Vt3=Vtl*sqrt (Tau) ; 

X3m=Xlm/ sqrt (Tau) ; 

Gammastar3=FNGammastar (X3m, Gamma) ; 

A3oal=  (Gammastarl /Gammastar3 )  *  (Kblock;/Kblock;3 )  *  (sqrt  (Tau)  )  /Prttavg; 
A3=A3oal*Al; 

g, 

o 

o 

o 


o 

o 

% - INPUT  FRACTIONAL  MEAN  RADIUS  CHANGE 

g, 

o 

Rfraction=input ( 'ENTER  FRACTIONAL  MEAN  RADIUS  CHANGE  FOR  NEXT  STAGE'); 
Rm3=Rm*Rf raction; 

H3=A3/ (2*pi*Rm3) ; 

Ho2rm=H3/ (2*Rm3) ; 

Rhub  tip3= ( l-Ho2rm) / ( l+Ho2rm) ; 

Um3=Um*Rf raction; 

Utip3=Um3* (l+Ho2rm)  ; 

Beta3m=atan (Um3/ (X3m*Vt3) ) ; 

Pt 3=Ptl* Prttavg ; 

Tt3=Ttl*Tau; 

PiO=Pt3/PtO; 

TauO=Tt3/TtO; 

Diatip=Rm+H/ 2 ; 

Diatip3=Rm3+H3/ 2 ; 

g, 

o 


Calculate  fan  performance  if  correctly  expanded  nozzle  was  added  here 


Xe=sqrt  (1-  (1/  (Pid*Pin*PiO)  )  (I/Gogl)  )  ; 
Thrust=Mdot*Xe*Vt3/ Gc; 

Hpower=Mdot*Cp*TtO* (TauO-1) *778/550; 
disp ( '  '  )  ; 

disp ( 'OVERALL  PRESSURE  RATIO  [T-T]  =  ' ) ; disp ( Pi 0 ) ; 

disp ( 'OVERALL  TEMPERATURE  RATIO  [T-T]=  ' ) ; disp ( TauO ) ; 
disp ( '  ' ) 

disp ('Fan  Diameter  = 

disp ('Fan  Static  Thrust  [=m*Ve/gc]  = 
disp ('Fan  drive  horse  power  = 

disp ( '  ' )  ; 

disp ( ' FOR  NEXT  STAGE : ' ) ; 
disp ('TIP  DIAMETER 
disp ( ' TIP  SPEED 
disp ('MEAN  RADIUS 
disp ( 'HUB-TIP  RATIO  = 
disp (' STAG. PRESSURE  = 
disp (' STAG. TEMP, 
disp ( 'BETA (MEAN-LINE) = 
disp  (  '  '  )  ; 

disp ( '  ' )  ; 

Xxx=input ( 'ENTER  I  FOR  ANOTHER  STAGE'); 

Utip=Utip3 ; 

Rhub  tip=Rhub  tip3; 

Ptl=Pt3; 

Ttl=Tt3; 


' ) ; disp (Diatip*2 ) ; disp ( '  inches ' ) ; 
' )  ; disp (Thrust) ; disp ( ' Ibsf ' ) ; 

' ) ; disp (Hpower ) ; disp ( ' HP ' ) ; 


)  ;disp (Diatip3*I2) ;disp ( '  inches ' ) ; 
) ;disp (Utip3) ;disp ( '  FT/SEC ' ) ; 

) ;disp (Rm3*I2) ;disp ( '  INCHES ' ) ; 

) ; disp (Rhub_tip3 ) ; 

); disp ( Pt3 ); disp ( '  PSIA'); 

) ;disp (Tt3) ;disp ( '  DEG.R'); 

); disp (Beta3m) ; disp ( '  DEG.'); 
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Be t lm=Be t a 3m ; 


Function  Calls : 

function  X=FNXfromm (M, Gamma) 

X= ( Gamma- 1) *M*M/2; 

X=X/ (1+X) ; 

X=sqrt (X) ; 

q, 

o 

function  Shock;ptratio=FNShock;ptratio  (Mwl ,  Gamma)  ; 

Msq=Mwl ^2 ; 

Gal= (Gamma+1 ) /2 ; 

Ga2= (Gamma-1 ) / 2  ; 

Ga3=2*Gamma/ (Gamma+1) ; 

Ga4=Gamma/ ( Gamma- 1) ; 

Ga5=l/ (1-Gamma) ; 

Shock;ptratio=  (  (Gal*Msq/  (l+Ga2*Msq)  )  "'Gai)  *  (  (Ga3*Msq-Ga2/Gal)  "'GaS) 

q, 

o 

function  Shock;pratio=FNShock;pratio  (Mwl ,  Gamma)  ; 

Msq=Mwl^2; 

Ga3=2*Gamma/ (Gamma+1) ; 

Ga6= (Gamma-1 ) / (Gamma+1 ) ; 

Shock;pratio=Ga3*Msq-Ga6 ; 

q, 

o 

function  Omegas tc=FNOmegastc ( Sig, B1 , B2 , Spch, Clrh) ; 

Binf=atan ( . 5*tan (B1 ) + . 5*tan (B2  )  )  ; 

Cl=2* (tan (Bl) -tan (B2) ) * (cos (Binf ) ) /Sig; 

Cdi=. 04*Cl*Cl*Sig*Spch+.25*Cl*Cl*Sig*Clrh/cos (B2) ; 
Omegastc=Cdi*Sig*  (  (cos  (Bl)  )  ''2)  /  (cos  (Binf)  )  '^3; 

q, 

o 

function  Omegap=FNOmegap (D, Sig,  Bl ,  B2  )  ; 

0megap=2 *Sig*  (  (cos  (Bl)  )  ^2)  *  ( .005+.16*D’^4)  /  (  (cos  (B2)  )  '^3)  ; 

q, 

o 

function  M=FNMf romx (X, Gamma) ; 

M=sqrt  (  (2/  (Gamma-1)  )  *X''2/  (1-X^2)  )  ; 

q, 

o 

function  Bet2=FNInv  dfact (D, Phi, Sig, Betl) ; 

A=(l-D+sin(Betl) / (2*Sig) ) /Phi; 

B=A^2+  (  (cos  (Betl)  )  '^2)  /  (4*Sig''2)  ; 

C= (1+sqrt (1- (1-1/ (4*Sig^2) ) *B/A^2) ) *A*cos (Betl) /B; 

Bet2=acos (C) ; 

q, 

o 

function  Gammastar=FNGammastar (X, Gamma) ; 

Gammastar=  (2*Gamma*X/  (Gamma-1)  )  *  (1-X''2)  (1/  (Gamma-1)  )  ; 

q, 

o 

function  D=FNDf act ( Phi , Sig, Betl ,  Bet2 )  ; 

D=l-Phi*cos (Betl) /cos (Bet2) + (tan (Betl) - 
Phi*tan(Bet2) ) *cos (Betl) / (2*Sig) ; 
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